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We address the question whether observables of an exactly solvable model of electrons coupled 
to (optical) phonons relax into large time stationary state values and investigate if the asymptotic 
expectation values can be computed using a stationary density matrix. Two initial nonequilibrium 
situations are considered. A sudden quench of the electron-phonon coupling, starting from the 
noninteracting canonical equilibrium at temperature T in the electron as well as in the phonon 
subsystems, leads to a rather simple dynamics. A richer time evolution emerges if the initial state is 
taken as the product of the phonon vacuum and the filled Fermi sea supplemented by a highly excited 
additional electron. Our model has a natural set of constants of motion, with as many elements 
as degrees of freedom. In accordance with earlier studies of such type of models we find that 
expectation values which become stationary can be described by the density matrix of a generalized 
Gibbs ensemble which differs from that of a canonical ensemble. For the model at hand it appears 
to be evident that the eigenmode occupancy operators should be used in the construction of the 
stationary density matrix. 
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I. INTRODUCTION 

It is of fundamental interest to reveal the conditions 
under which an isolated quantum system for t cx) re- 
laxes into a state that can be described by a stationary 
density matrix. Furthermore, a detailed understanding of 
the relaxation process is desirable. We here characterize 
relaxation by considering the time evolution of expecta- 
tion values of observables. Starting out from a nonequi- 
librium situation at t = described by the initial density 
matrix p{ we ask whether the expectation value 



Tr p(t)Oi 



(1) 



of a local observable 0\ for time t — )> oo approaches a 
constant value which can also be computed considering 
the stationary density matrix p^t 



lim (O 



Tr 



PstOl 
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Here p{t) denotes the statistical operator at time t which 
follows from solving the von Neumann equation for the 
given initial condition p{t = 0) = pi. A "local observ- 
able" is defined as one which only contains degrees of 
freedom from a subsystem S of the isolated quantum sys- 
tem Q. We focus on measurements in subsystems and one 
must thus be careful in interpreting pst as the density ma- 
trix describing the stationary state of the entire quantum 
system Q. For a meaningful description of the stationary 
state Pst should be independent of the chosen observable. 
To avoid recurrence effects one has to perform the ther- 
modynamic limit of Q which is often done by taking the 
thermodynamic limit Vf — >■ oo, with the volume Vs of 
the environment £ = Q/S, keeping Vg fixed. 

Relaxation to a time independent expectation value in 
the strict sense can only occur after the thermodynamic 



limit has been taken. Alternatively one can address the 
question whether a "quasi stationary state" is reached 
in a finite system. By this we understand a situation in 
which expectation values "fluctuate" around a constant 
value which can be extracted by averaging over time. If 
such a state is reached one can ask if the time averaged 
value can be computed using a stationary statistical en- 
semble. In the averaging it might even be meaningful to 
increase the time interval beyond the characteristic time 

= L/v — denoted recurrence time in what follows — 
with L being a typical length and v a typical velocity of 
the system (see below). 

In equilibrium statistical physics we commonly work 
with thermal ensembles as the ones describing the state. 
They are characterized by the density matrix 

p= - e- ^7=1 ''^■fj- , Z = Tre- ^?=i "^^'^ , (3) 

with the partition function Z({rij}). The sum usually 
runs over only a few terms containing operators Ij , such 
as the Hamiltonian H and the particle number operator 
N corresponding to the macroscopic variables energy and 
particle number. The Lagrange multipliers rjj are fixed 
such that the expectation values of the Ij take given val- 
ues (e.g. given average energy and particle number) 



r(0) 



Chosing the corresponding Lagrange parameters in Eq. 
([3]) maximizes the entropy S = Tr[/51n/5~^] under the 



constraint of fixed /, 



(0) 1 



Within the observables of a 



closed system (with fixed particle number) the Hamil- 
tonian plays a special role as in many situations the ex- 
pectation value of the energy is the only conserved quan- 
tity. The corresponding Lagrange multiplier is the in- 
verse temperature /? = 1/T and the thermal ensemble 



2 



with only the energy expectation value fixed is the canon- 
ical one. 

Jaynes- studied generalized ensembles — now com- 
monly referred to as generalized Gibbs ensembles 
(GGEs) — in which additional observables Ij besides the 
energy are assumed to take a given expectation value 
(fixed if''^) and together with the corresponding La- 
grange multipliers enter the sum in Eq. ([3]). 

Starting out with the initial nonequilibrium state given 

by Pi and under the assumption that {0\) converges 

\ / p(t) 

for t — > oo one might expect that the stationary expecta- 
tion value can be computed from Eq. ([2]) with pst = Pcan 
and 



Pc 



(4) 



This is known as thermalization. The inverse tempera- 
ture /3 is set by the constraint 



= H 



H 



In particular, this is the expected behavior if the energy 
is the only (independent) constant of motion. 

Recent experiments in the field of ultracold atoms^ led 
to a revived interest into relaxation dynamics. Due to 
the long coherence time such systems are ideal candi- 
dates to study the relaxation into a stationary state in 
a controlled setup. In the experiments an equilibrium 
state is disturbed by a sudden quench of system param- 
eters. The experiments led to several interesting theo- 
retical studies in which models, usually considered in the 
field of quantum many-particle physics, were investigated 
concerning their relaxation properties. In analogy to the 
experiments in most studies the system is assumed to be 
in an eigenstate (e.g. the noninteracting groundstate) or 
a canonical thermal state with fixed temperature T of 
the initial many-body Hamiltonian Hi. At t — Q model 
parameters, in most cases the two-particle interaction be- 
ing a crucial element of the models, is quenched instan- 
taneously to a different value and the dynamics of the 
initial state under the time evolution given by the final 
Hamiltonian Hf is computed. 

The investigations can be grouped in three classes. 
Analytical studies of models which can be solved 
exactlyr'"— numerical studieSfi^"— and approximate an- 
alytical studiesii^iiS For most of the considered models 
the relaxation properties of a restricted set of observables 
was computed. A more general perspective on theprob- 
lem for a certain class of models is taken in Refs. [20 and 
,21. using methods of boundary critical phenomena and 
conformal field theory. The ultimate goal of the "case 
studies" is (i) to derive criteria which a priori allow to 
answer the question whether certain (local) observables 
become stationary and (ii) to construct the stationary 
density matrix pst, which is the appropriate ensemble 
by which the asymptotic expectation values can be com- 
puted. It is of particular interest to understand the con- 
ditions under which the latter becomes the canonical one 



and the system thermalizes. From our considerations it 
is plausible to expect that the number and character of 
the constants of motion of a specific model are of crucial 
importance in answering the questions addressed above. 
Roughly speaking, if the number of constants of motion 
is large the time evolved state contains a lot of infor- 
mation from the initial state and thermalization cannot 
be expected. In all the studies mentioned above single 
component systems containing either bosons or fermions 
were studied. Besides the asymptotic long time behav- 
ior transient nonequilibrium effects such as "collapse and 
revival" were investigated. 

From the analytical^i^iii2r2ii2S and numerical— il stud- 
ies of exactly solvable models increasing evidence was col- 
lected that if the expectation value of an observable O 
approaches a constant large time limit, the latter cannot 
be obtained using the thermal canonical density matrix. 
Instead stationary density matrices of the GGE type Eq. 
([3]), corresponding to situations with more restrictions 
(than a fixed average energy) set by the initial state hav- 
ing a density matrix 



1 



Pgge 



Zgge 



=- ^.7=1 v,i, ^ Zgge = Tr e' ^^=i "^^^ (5) 



turned out to be promising candidates. On the basis 
of the above considerations this is not surprising, as the 
models considered are characterized by more integrals of 
motion than just the energy and one thus expects for 
the stationary state an increased "memory" of the initial 
state. A set of constants of motion were taken as the 
Ij and the Lagrange multipliers rjj were determined such 
that 



PGGE 



(6) 



The choice of conserved observables is not unique. E.g. 
with H being a constant of motion the same holds for 
and it was shown^^ that under certain conditions 
the GGE expectation value of a given observable O 



O) = Tr pggeO 

PGGE 



(7) 



might depend on the selected set. Thus the GGE de- 
scribes the stationary state only if the "correct" set of 
"independent" Ij is chosen. For a general model with 
several possible sets of constants of motion up to now it 
is not clear how to select the "correct" set a priori. 

In cases in which the system's Hamiltonian can be 
brought into the form 



(8) 



with either fermionic or bosonic creation and annihilation 
operators '0^^'' and quantum number k the eigenmode oc- 
cupation operators fik = constitute a natural set of 
conserved quantities, with as many elements as degrees 
of freedom. Then the statistical expectation value taken 
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with the PGGE set up with Ik = fik and Lagrange mul- 
tiphers fixed by the constraints Eq. © led for generic 
model parameters to the correct large time limit of (cer- 
tain) observables provided the latter exists^iSiiii and the 
initial states are homogeneous.^^ 

Cold atom gases is not the only subfield of modern 
condensed matter physics in which understanding the re- 
laxation dynamics is of crucial importance and quenches 
are not the only process leading to an interesting time 
evolution. In the area of photoexcited semiconductors 
much effort has been put into measuring (pump-probe 
techniques) and understanding the relaxation dynamics 
of highly excited electrons coupled to optical phonons; for 
a recent review see Ref. .23. . The long-time asymptotics 
of the electron-phonon system is difficult to study ex- 
perimentally due to the strong coupling to other degrees 
of freedom (e.g. acoustic phonons and holes), but the 
short-time (transient) dynamics shows interesting non- 
Markovian effects requiring a treatment beyond the use 
of Boltzmann equations<^"— A detailed understanding of 
relaxation also plays a major role in the field of condensed 
matter system based quantum information processing?21 
Here the strong coupling of the degrees of freedom envis- 
aged as quantum-bits to the environment usually leads 
to a coherence time to short to perform substantial infor- 
mation processing. Gaining insights into the relaxation 
process might lead to ways to circumvent this obstacle 
and significantly increase the coherence time. 

We here supplement the recent "case studies" on the 
relaxation dynamics of either bosonic or fermionic cor- 
related systems by analytically investigating the time 
evolution of a two- component model of electrons cou- 
pled to phonons. The model naturally contains two 
subsystems — the electron and the phonon systems — and 
an observable can be considered as local if it contains 
only fermionic or bosonic degrees of freedom. We first 
consider the time evolution resulting from a quench of 
the electron-phonon coupling from zero to a finite value 
starting with the noninteracting canonical equilibrium at 
temperature T in the electron as well as the phonon sub- 
systems. In addition we study the dynamics inferred by 
the interacting Hamiltonian out of a pure state given by 
the product of the phonon vacuum and the filled Fermi 
sea supplemented by a highly excited additional electron 
of momentum fcg . The short-time dynamics of our model 
starting with this initial state was earlier discussed in the 
context of optically excited semiconductors^^ and used 
to explain results of pump-probe experiments i^i^^ While 
the time evolution of observables, in particular the elec- 
tron or phonon momentum distribution function and the 
subsystem energies, is rather simple for the quench a rich 
dynamics is found in case of the "fco-excitation" . We 
show that in the large time limit t — >■ oo the subsystem 
energies (and the energy in the electron-phonon coupling) 
converge to stationary values for both nonequilibrium ini- 
tial states. The same holds for the electron momentum 
distribution, while for the phonon momentum distribu- 
tion function convergence is only achieved after averag- 



ing over a small momentum interval. Our model can be 
brought into the form Eq. ([H]). It thus contains (at least) 
as many constants of motion as degrees of freedom. In 
accordance with earlier studiea^i^iii we find that for both 
initial states the expectation values of observables which 
become stationary can be described by the density ma- 
trix of a generalized Gibbs ensemble with the eigenmode 
occupation operators chosen as the Ij. 

The rest of this paper is organized as follows. In 
Sect. |TT]we introduce our one-dimensional (Id) electron- 
phonon model and specify our initial states. The 
bosonization of the fermionic field operator, which al- 
lows to obtain analytical results for the time evolution 
of fermionic observables, is discussed in Sect. IIIIl Sec- 
tion IIVI is devoted to the GGE for the model at hand. 
In Sects. |V] and IVII we discuss the relaxation dynamics 
for our two distinct initial states and compare the long- 
time expectation values to those obtained from the GGE. 
Finally, our results are summarized in Sect. IVIIl 

II. THE ELECTRON-PHONON MODEL AND 
THE NONEQUILIBRIUM INITIAL STATES 

We consider a model of electrons on a Id ring of length 
L (periodic boundary conditions) coupled to phonons by 
a Holstcin type electron-phonon interaction given by the 
Hamiltonian 

H = H, + Hp + H,p, (9) 

with 

He = E^'^(4«fe-(4a/.)J . (10) 

k 

Hp = J2^,BlB^, (11) 

q>a 

^-P = [tJ E5(9)(4+A^, + H.c.) .(12) 

q>0 

Here cik (Bq) is the annihilation operator of an electron 
(phonon) with momentum k (q) and (...)g is the expec- 
tation value in the noninteracting groundstate (normal 
ordering). The momentum dependent electron-phonon 
coupling is denoted by g{q). To be specific we consider 
the (simple) form g{q) = g&{qc — q)- We restrict our 
considerations to a single branch of chiral (right-moving) 
spinless fermions with linear single-particle dispersion 
Efc = VF{k — k-p), the Fermi velocity v-p, and the Fermi 
momentum fcp which, without loss of generality, is set 
to zero in the following (fcp = 0). We assume that the 
fermion states do not have a lower bound (Dirac model) 
and all momentum states with fc < are filled in the 
groundstate (filled Fermi sea). For the sake of conve- 
nience the Fermi momentum fcp corresponds to the first 
empty instead of the last occupied state. Divergencies 
possibly resulting from these states are regularized by 
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the normal ordering. Furthermore, we focus on optical 
phonons (Einstein model) with a single energy ujg — ujq. 

The assumed g-dependences of the phonon dispersion 
and the electron-phonon coupling can be relaxed with- 
out spoiling the possibility of an exact analytical solu- 
tion of the model. Varying the function g{q) only leads 
to minor changes in the short-time dynamics^^ (as long 
as g{0) remains finite) and we do not see any physical 
reason why our main results for the long-time relaxation 
obtained here should dependent on the precise form of 
g{q). Our model is a variant of the purely fermionic 
Tomonaga-Luttinger mode P^i^^ — more precisely what is 
called the chiral g4-model2^ — and similar to this case the 
linear fermion dispersion is crucial for an exact analytic 
solution. Equilibrium properties of the above model, 
also considering acoustic instead of optical phonons, 
were discussed earlier.—"— The quench dynamics of the 
Tomonaga-Luttinger model was studied in Ref. '5'. 

To diagonalize our Hamiltonian Eqs. (P))- ([T^ we first 
introduce the electron density operator 



aia 



k"'k+q J 



(13) 



with g > 0. For q = we define the electron number 
operator relative to the groundstate (filled Fermi sea) as 



N = 



k 



(4«fc-(4«fe)J ■ (14) 



With a proper normalization the densities dq obey Bose 
commutation relations. If one defines 



for (7 > 0, the commutation relations read 



1/2 



= 6, 



9.9' 



= 



(15) 



(16) 



The electron-hole excitations in H^p can straightfor- 
wardly be written linearly in the boson operators and 
&|. For Id systems with a linear dispersion it is in addi- 
tion possible to write the kinetic energy of the fermions 



„29.36 



g>0 



(17) 



where c(iV) only contains the fermionic particle number 
operator. In the following this term can be droped be- 
cause we are working in a sector of the Hilbert space 
with constant particle number. We can now formulate 
Eqs. (P|)- p^ in terms of the bosonic fermion densities 
and the phonons as 



H 



^VFqblbg 

+9 



q>0 



H.c 



(18) 



From this expression it becomes evident that the ground- 
state of the model is still the tensor product of the 
Fermi sea — corresponding to the vacuum of the bosonic 
fermion density bq — and the phonon vacuum. Therefore 
all groundstate expectation values, like e.g. the fermionic 
momentum distribution, are given by the noninteracting 
ones. 

Using a canonical transformation the problem of the 
coupled bosonic modes Eq. (IT51) can be brought in the 
form of Eq. ([Sl).'^^ The transformation is given by 



B„ 



^qCq (^q^q 
^q^q ^ Pq^q 



with 



2 ^ |A+(g) -(JqI 2 

^ A+(q)-A_(g) ' 

and the mode energies 



|A-(g) - ^ol 
A+(g)-A_(g) 



(19) 



(20) 



A±(<?) = ^ |wFg + wo± \/ {vpq - ujof + Ag'^qQiqc ^ '?)| 

(21) 

Note that c^-f = 1 for all 9 > and that = 1, ^ 
for q > qc- Here we focus on the case wf9c > f^o- In the 
new bosonic operators ctq and /3g the Hamiltonian reads 



H = 



^-{qWq 



(22) 



(?>0 



0<g<gc 



In order to obtain a stable groundstate the boson energies 
\±{q) have to be larger than zero4i For a given wq this 
leads to a restriction of the coupling strength g that can 

be used. The dimensionless parameters of the model are 

2 

r = — — for the electron-phonon coupling, VI — 

v-pqc ^ 05 v-pqc 

for the phonon frequency, and v = for the inverse 
of the ring length. Stability requires that F < 17. The 
momentum dependence of the eigenmode energies A± (g) 
and the coefficients and are shown in Fig. [1] for a 
typical set of parameters with F = 0.01 and = 0.1. 

With Eqs. ((22|) and ((20| computing the dynamics of the 
phonon ladder operators and fermionic densities becomes 
simple. We find 

Bqit) = CqSq (e-'^+(')* - g-^^-^')*) bq 

+ (s^e-'^+(«)* + c2e-'^-(«)*)s, (23) 

and 

bqit) = (c2e-^^+(^)* + s2e-^-(^)*)6, 

+CqSq (e-'^+(9)* - e-'^-(«)*) Bq . (24) 

The time dependence of expectation values of observables 
which can be written in terms of the Bq^'^ and fog^'' can 
thus be expressed by expectation values taken with the 
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FIG. 1: (Color online) Eigenmodes A±(g) and coefficients 
and of the eigenvectors for F = 0.01 and Q. = 0.1. Note the 
sharp transition (discontinuity) to the noninteracting values 
at Qc. 



initial density matrix p\ . This gives us direct access to the 
dynamics of the phonon momentum distribution function 



think of this state being (approximately) realized in a 
doped semiconductor (Fermi sea in conduction band) in 
which an additional "hot" electron with momentum fco is 
optically pumped from the valence band into the conduc- 
tion band. Both initial conditions correspond to nonequi- 
librium states if the time evolution is performed with the 
Hamiltonian Eq. ([TSl) for 5^0. We note in passing that 
starting out with the noninteracting groundstate (T = 0) 
would not lead to a time dependence of expectation val- 
ues after a sudden quench of g as in the present model 
this state remains the groundstate even for g ^ 0. 



III. BOSONIZATION OF THE FERMIONS 

To calculate fermionic expectation values like the mo- 
mentum distribution function we use the bosonization of 
the fermionic field operator 



(29) 



(25) One can prove the operator identitys^ 



and that of the subsystem energies for both nonequilib- 
riuni states considered. Alternatively with Eqs. (1^51) and 
(1241) the time evolution of the density matrix can be given 

in a closed form if p\ can be expressed in terms of the feq^"* 
and B'^q^ (see below). To compute the fermionic momen- 
tum distribution function 



n(fc, t) 



al.(0)afc(0) 



m 



one has to establish a relation between the fermionic op- 
erators ajj^'' and the h^q^ . This bosonization (of the field 
operator) will be discussed in Sect. IIIIl 

We consider two different initial situations at time 
i = 0. In the first one the decoupled {g = 0) electron- 
phonon system is initially assumed to be in a thermal 
state with a common temperature T = of the elec- 
tron and phonon subsystems corresponding to the non- 
interacting canonical ensemble. It is determined by the 
initial density matrix (superscript q for "quench" ) 



with 



— ixTr/L , 



(30) 



(31) 



where TJ^ denotes a unitary fermionic raising operator 
which commutes with the fog^'' and maps the A^-electron 
groundstate to the {N + l)-electron groundstate. As we 
are interested in the dynamics for a fixed particle number 
neither the term proportional to N in Eq. (1311) nor the 
fermionic raising operator W affect the result and both 
can be dropped in the following. Using Eqs. (|3T|). ((30|) 
and the back transform 



1 



L/2 



L/2 



e''=^^/'^(a;) dx 



/z 



(26) of Eq. (1^ as weU as Eq. ((M| the time dependence of a 



The second initial density matrix is given by a pure state 
(superscript fco for fco-cxcitation) 



(27) 



can be expressed in terms of the and Bq'. As further 
deepened below this allows us to give a closed expression 
for the time evolution of the fermionic momentum dis- 
tribution n{k,t) as well as its value in the appropriate 
GGE. 



?(t) 



(t) 

k 



with 



i*i) = «IoiFs; 



vac) 



(28) 



IV. THE GENERALIZED GIBBS ENSEMBLE 



Here |FS) denotes the filled Fermi sea (vacuum with re- 
spect to the bq) and |vac) the phonon vacuum. One can 



As shown in Sect. |TT] our coupled electron-phonon 
Hamiltonian can be brought into the diagonal form Eq. 
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([8]) with the bosonic operators OL^q ^ and /3q^'' and the cor- 
responding mode energies \±{q). Therefore the eigen- 
mode occupancies djdq and /3j/3g constitute a set of con- 
stants of motion which has as many elements as degrees 
of freedom in our modeL The occupancies thus form a 
natural set of operators Ij which can be used to set up 
the density matrix of a GGE Eq. ([5]). We again empha- 
size that in Ref. [Tl| an example is given, which shows that 
an alternative choice of conserved Ij might lead to GGE 
expectation values of observables which differ from the 
ones obtained by the natural choice. We show here that 
for our model all the studied observables which become 
stationary in the long time limit approach a value con- 
sistent with the GGE set up by the set of occupancies 
(called natural GGE below). 

The GGE is described by the density matrix 



Pgge 



Zgge 



(32) 



with the Lagrange multipliers "qq and ^q determined by 
the initial condition 



Pgge 



(33) 



and similarly for j3q. For a density matrix of the form 
Eq. (15^ the eigenmode occupancies are given by^ 



{a\aq) = n^{T^q) , (/3t/3,\_ ^ n^{iq) (34) 

^'GGE \ 1 / pQCE 

with the Bose function 

To fix the Lagrange multipliers we still have to com- 
pute the left hand side of Eq. ([55]) for the two initial 
density matrices Eqs. (f26|) and ([27l) . For the quench we 
obtain using (the inversion of) Eq. ([Tl 



'q^i/ a ^ s nB(^WFg) +c nB(/3wo) , (35) 
Pi 



which leads to the set of nonlinear equations for the rjq 
and £,q 

nBiCq) = s^nB(/3wF'?) + c^?^B(/3l^o) • (36) 

Obviously, the rjq and ^q are functions of the inverse tem- 
perature (3. 

For the fco-excitation it follows using (the inversion of) 
Eq. (jl9p and the "vacuum" -properties of the initial pure 



state that 



= Tr 



Pi'alaq 



(vac| (g) (FS I auo (cg^J 

X {cqhq + SqB^ fl^^ |FS ) (g) |vac) 

cl {YS\ak,b\bqal |FS) 



27r 
Tq 
2tt 
Tq 



(FS|afc„_,4^_^|FS 



e(fco - q) , 



(37) 



with the 0-function defined such that 9(0) = 1. In the 
step from the fourth to the fifth equation we have used 
twice that from Eqs. ([T^ and it follows that 



'■ko-q 



Similarly we obtain 



2tt 



(38) 



which determines the Lagrange parameters using Eqs. 
(131 and (El. 



A. Expectation values for the quench 

We are now in a position to determine the GGE expec- 
tation values for the subsystem energies, the phonon mo- 
mentum distribution NQQ^{q), and the fermion momen- 
tum distribution nGGE(fc) first focusing on the quench. 
With Eqs. (0, and ^ we straightforwardly ob- 

tain 



PGGE 



^ vpq {uBiPvpg) - 2Cg s 
g>0 

X [uBiPvpg) - UBi/Sujo)]} 
Wo ^ {uBiPuo) + 2c^ si 

g>0 

X [uBiPv-pq) - 7iB(/3a;o)]} 

25X1 V^CgSg (c^ - si) 
g>0 

X [nB{l3vFq) - nB{/3ujo)] ■ 



(39) 



The first terms in the subsystem energies are the expec- 
tation values taken with the initial canonical ensemble at 
g = 0- Without a high momentum cutoff the energy in 
the phonon subsystem obviously diverges. Here we avoid 
to introduce such a cutoff by considering the excess ener- 



gies 6 (Hp 



PGGE 



Hr 



PGGE 



Hr 



and S (He , 



PGGE 
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(defined similarly) in the following. Note that the mo- 
mentum sums containing a factor Sq are cut off at as 



for g > gc- With 



q>0 



dq . 



for large L it becomes apparent that the subsystem excess 
energies and the energy in the electron-phonon coupling 
are extensive and scale ~ L. 

For the phonon momentum distribution function it fol- 
lows similarly that 



riBiPvFq) - UBifioJo)] . (40) 



It is instructive to compare this result to the canonical 
(equilibrium) phonon distribution of the interacting sys- 
tem (same parameters g and ujq) at some temperature 
T — 1/(3 characterized by the density matrix Eq. A 
straightforward calculation using the methods introduced 
above gives 

^can(g) - slnB0X+iq))+clnB0X-iq)) • 

Obviously, the g-dependences of this function and Eq. 

differ and no temperature T can be found leading to 
coinciding results. 

Using the method introduced in Sect. Illli the Baker- 
HausdorfF relation, and the formulai^ 



gAgB\ ^ ^{A-'+2AB+B'}/2 



one obtains for the fermionic momentum distribution 
function in the GGE 



»^GGE(fc) 



27r 



-iqix — iO) 



(41) 



1 /"^/^ 

- dxe exp < > 

X exp ^ [1 - cosiqx)] UBiPypq)^ 

X exp 1^ ^ cpl [1 - cos(ga;)] [71b(^wf9) - nB(/3a;o)]| 

The third factor, which vanishes for g = 0, contains the 
information about the electron-phonon coupling. The 
first two terms constitute the canonical momentum dis- 
tribution at temperature T = for a noninteracting 
fermionic system with linear dispersion on a ring of size 
L in equilibrium. This g = expression was earlier de- 
rived in Ref. [H. As discussed there it becomes equal 



to the (grand canonical) Fermi function [e 



1]" 



(only) in the thermodynamic limit L — )■ cx). Instead of 
numerically performing the sums and the integral for fi- 
nite L Eq. (HIT) can very efficiently be evaluated using an 
iterative approach introduced in the Appendix of Ref. Issl . 



Adopting this method to the present situation we obtain 
for m S Z, km = 27rm/L, I — 0,1,2, . . . 



"-GGE(fcm) = ^ fln+m , 



(42) 




1=1 



f{l) = nBiPvpqi) - 2cl^sl^ [nBiPvFqi) - nB(^Wo)] 



with qi = 27:1/ L. As for the phonons we compare this re- 
sult to the thermal distribution with temperature T. The 
canonical fermion distribution function of the interacting 
system can be computed along the same lines as the GGE 
distribution. Because of the involved structure a compar- 
ison of the analytical expressions is less instructive as for 
the phonons. In Fig. [5] we therefore compare numerical 
results for the GGE with T = 0.01, VI — 0.1, system size 
parameter v — 2tt /{Lqc) — 10"'^, and dimensionless tem- 
perature r — T /{v-pqc) = 0.1 with the best fit of a canon- 
ical distribution function for the same parameters — in 
particular the same electron-phonon coupling and the 
same system size — and fitting parameter f = T /(v-pqc) . 
The best agreement is achieved for fb = 0.10275. In gen- 
eral Tb depends on the model parameters and t. The 
differences are small but significant as becomes explicit 
in the inset which shows the absolute value of the differ- 
ence of the two distributions. Doubling the system size 
(halving v) does not lead to any changes on the scale of 
the main plot as well as the one of the inset. Thus the 
curves can be considered to be in the thermodynamic 
limit and the two ensembles lead to different results even 
after the latter has been performed. This is a crucial ob- 
servation as we have to distinguish this type of deviation 
between the predictions of two ensembles from the one 
which might appear at finite L but vanishes for L ^ oo. 
An example for the latter case is the difference between 
the canonical and the grand canonical ensembles as re- 
ferred to in lectures on statistical mechanics. For non- 
interacting fermions with a linear dispersion such finite 
size differences are explicitly studied in Ref. Issl . 



B. Expectation values for the fco-excitation 



We next derive the same expectation values but for the 
/co-excitation. Using Eqs. (HH), dHT]), and ^ we 
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FIG. 2: (Color online) The GGE and canonical momentum 
distribution function of the fermions. The system parameters 
are F = 0.01, Q = 0.1, u = 27r/(Lgc) = 10"^ and the di- 
mensionless initial temperature is r = T/{vFqc) = 0.1. The 
canonical distribution is the best fit to the GGE one with 
the fitted temperature fb = 0.10275. The inset shows the 
absolute value of the difference of the two functions. 



obtain 



PGGE 



PGGE 



PGGE 



X! — (cg + 4) ©(fco - q) , 

i3>0 ^ 



(43) 



25X1 ]r V^CqSq [cl - si) e(fco - q) . 

q>0 ^ 



In contrast to the quench case the energies are not exten- 
sive. This is related to the fact that even for L — >■ 00 we 
only add a single additional fermion at momentum to 
the filled Fermi sea. The phonon momentum distribution 
function is given by 



27r 

iVGGE(9) = 2— slQ{ko-q) 
Lq 



(44) 



Using the bosonization of the fermionic fields the GGE 
expectation value for the momentum distribution func- 
tion follows as 



1 r^l'^ { 27r 

"GGE(fc) = Y / dxe-*'^'' exp <^ V — 



-iq{x — iQ) 



X exp 1 2 ^ (^^^ (4 + 4) [1 - cos(ga;)] e(fco - -?) | • 

It is now crucial to realize that due to the factor 1 /L^ in 
the exponent the second term approaches 1 in the ther- 
modynamic limit L — >■ 00 independent of the electron- 
phonon coupling. In this limit the remaining terms form 
a step function. We thus find 



lim nGGE(fc) = B(fc) 



(45) 



and the fermionic momentum distribution function of the 
GGE in the thermodynamic limit becomes equal to the 
one of the groundstate (which for the present model is 
equal to the noninteracting one; see above). This is con- 
sistent with the observation that the energies for the k^- 
excitation are not extensive. For finite L, nQQ^[k) can 
again be computed iteratively using Eq. (j42|) with f{l) 
replaced by 



Q{ko - qi)/l 



(46) 



We note in passing that as for the quench for finite 
L, nQQ-E{k) is different from the canonical distribution 
function obtained for the same system parameters at an 
optimally chosen temperature T. The same holds for 

^^GGE(g)- 



V. THE QUENCH DYNAMICS 

We now investigate the dynamics of the density matrix 
and the expectation values considered above under the 
Hamiltonian Eq. ([9]) starting out from the noninteracting 
{g = 0) canonical density matrix Eq. (|26|) at temperature 
T — 1//3. Using the formal solution 



(47) 



of the von Neumann equation and Eq. (|19|) we end up 
with 



Zp'i(t)=exp <~f3^^{vFqcl+ujQsl)alaq 

[ q>0 

+c,Sg (wo - t'pg) e-'^^^^^'dt^, + H.c. 



{vpqsl + ujqcD I3l/3g 



(48) 



where 



AX{q) = X+{q) - X^{q) . 

We emphasize that only the difference of the two eigen- 
mode energies enter the dynamics. To compute the ex- 
pectation values of interest with p^{t) one can now diag- 
onalize the time dependent 2x2 coefficient matrix of the 
quadratic form appearing in the exponent of Eq. ()48|) and 
introduce new time dependent bosonic operators as lin- 
ear combinations of the dq^'' and /3q^\ Alternatively, one 
can use Eqs. ([23|) and ([24|) and compute the expectation 
values with p?. 

For the time dependence of the subsystem excess en- 
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ergies and interaction energy this leads to 



S(H,) = -yvFq2 



2 2 



X [riBiPvYq) - nB(/3wo)] [1 - cos{AA(g)t}] , 



, 2 2 
'-Cq Sq 



q>Q 

X [nBil3vFq) - nB(/3a;o)] [1 - cos{AA(g)i}] , 



CqSq (Cq Sg) 



g>0 



X [riBiPvYq) - nB(/3wo)] [1 - cos{AA(g)t}] .(49) 

It is easy to show that for all t the total excess energy 

5(H) sums up to zero due to energy conservation. 

/p''(t) 

To answer the question if the (excess) energies become 
stationary in the long-time limit we first perform the 
thermodynamic limit. Afterwards the oscillatory terms 
average out for t — >■ oo when the momentum integrals are 
performed. For L — > oo the (excess) energy expectation 
values (per particle) thus become stationary and equal to 
the GGE expectation values (per particle) as determined 
in Eq. ([39]) . This provides a first indication that for our 
model long-time expectation values of observables which 
become stationary can indeed be computed using the ap- 
propriate GGE. As discussed in the introduction in per- 
forming the thermodynamic limit one often keeps the size 
of one of the subsystems (the S) containing the relevant 
local observables fixed while the thermodynamic limit is 
performed in its complement, the environment E^ln our 
model the limit L — > oo is simultaneously performed in 
the fermion and phonon subsystems (quantization of mo- 
menta) . 

As shown in Fig. [3] (solid lines) for large t the energies 
oscillate around their asymptotic GGE values (dashed 
lines) with seemingly a single frequency and a slowly de- 
caying amplitude. From the data it is obvious that a 
rather good approximation for the stationary expecta- 
tion values can be obtained by averaging the (excess) en- 
ergies over a time interval which is much larger than the 
inverse oscillation frequency and starts at a sufhciently 
large time. 

Analytic insights on the long-time behavior can 
be gained by applying the techniques of asymptotic 
analysis such as the stationary phase method to the 
momentum integrals in Eq. Using the latter one 

can show that for not too strong couplings 2V < Q. (on 
which we mainly focus) such that the stationary point 
(of AA((7)) lies inside the integration interval (0,qc) the 
dominant oscillation frequency (at large times) is given 
by the minimum of the mode energy difference 



AA(gn 



2r 



1 



and the amplitude falls off as l/^/t."^^ The convergence 
towards the stationary values is thus rather slow. Fur- 
thermore, the second (high) frequency visible at small 
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FIG. 3: (Color online) Time evolution of the (excess) energies 
in the fermion (solid line, top) and phonon (solid line, middle) 
subsystems as well as of the energy in the interacting part 
of the Hamiltonian (solid line, bottom). The dashed lines 
for t £ [1000, 1500] are the corresponding GGE expectation 
values. The parameters are V = 0.01, 57 = 0.1 and the initial 
temperature is given by r = 0.1. The symbols in the inset 
show data for the energy in the fermion subsystem obtained 
at a finite system size v = 2 x 10~^. The recurrence time is 
then given by {vpqc)tr — vy2ii /(Lv)L/vy ~ 2-k/i' ~ 3142. 



times results from a contribution of the boundaries of 
the momentum integrals and can be identified as AA(gc)- 
The mode energy difference of the lower boundary does 
not appear as a frequency as the corresponding ampli- 
tude vanishes. 

For comparison the inset of Fig. [3] also contains results 
for the energy in the electron subsystem obtained at finite 
L (symbols). Obviously on the scale t^ = L/v-p finite size 
effects set in and the oscillations become rather erratic. 
Remarkably, for upper times not too large compared to 

the L —> oo asymptotic value can still be extracted 
accurately by averaging over the time even beyond t,-. 
The same holds for the other energy expectation values. 

From the energy in the phonon subsystem Eq. p9)) 
the time evolution of the phonon momentum distribution 
function out of the initial Bose function nB(/3i^o) can be 
read off 

6Nq{q,t) = Nq{q,t)^nBil3LUa) 

X [l-cos{AA(g)i}] • (50) 

As the phononic annihilation and creation operators are 
linear combinations of the eigenmode ladder operators 
(see Eq. (jl9|) ) it is not surprising, that the phonon mo- 
mentum distribution function for a fixed q does not be- 
come stationary; it shows a sinusoidal oscillation with 
frequency AA(g) and fixed amplitude for all t. The time- 
dependent part SNq{q, t) is shown in Fig. |3]for F — 0.01, 
n = 0.2, and temperature r — 0.1. The shape for fixed 
t can be understood from Fig. [T] (in which and are 
shown) and the g-dependence of the difference of the two 
Bose functions. Only after averaging over a small interval 
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FIG. 5: (Color online) Momentum distribution function of 
the fermions as a function of k for different t. The GGE 
distribution obtained for the same system size is shown for 
comparison. The parameters are F = 0.03, Q, — 0.1, ly = 10""^, 
and the initial temperature is given by r = 0.1. 



FIG. 4: (Color online) The time-dependent part of the phonon 
momentum distribution function 5Nq{q,t) as a function of q 
and t. The parameters are F — 0.01, — 0.2, and the initial 
temperature is given by r = 0.1. 



Aq around q, Nq{q,t) becomes stationary at large times 
as the cosine term drops out. 

This behavior has to be contrasted to the one of the 
fermion momentum distribution function nq(fc, t). By the 
bosonization of the field operator (see Sect. IIIip the rela- 
tion between the a^}"' and the eigenmode ladder operators 
is highly nonlinear and it is not clear a priori whether or 
not nq(k,t) becomes stationary. Using the methods in- 
troduced above we obtain 



nq(fc,t)==y / dxe ''"^expiV'^ 
X exp S - T- [1 ~ cos(gx)] UBiPvpg) > 
X exp i ^ 



-iq{x—iO) 



(51) 



[1 - cos{qx)] [uBif^vpg) - nB(/3wo)] 



X [1 - cos{AA(g)t}] 



First taking the thermodynamic limit and subsequently 
the long-time limit the oscillatory time dependent term 
again averages out and we find convergence to the GGE 
expectation value Eq. (|4T|) . 

lim lim nq{k,t) — lim nGGE(fc) ■ 

Furthermore, an analytic stationary phase analysis simi- 
lar to the one performed for the energies shows that for 
sufficiently large t and fixed fc, nq{k,t) oscillates with 



frequency AX{q^in) around nQQ^ik) with an amplitude 
which decays as 1/Vt. As for the energies the second 
relevant frequency is AA((jc). 

For finite (but large) L, nq{k,t) can again efficiently 
be computed numerically using the recursion relation Eq. 
(H^ with f{l) replaced by 



/(/) = UBiPvFqi) - 2c s [nBif3vFqi) - nBil^uJo)] 



X [l-cos{AA(gi)i}] 



(52) 



In Fig. [5] the time evolution of the fermion momentum 
distribution function as a function of k and a few t is 
compared to the GGE prediction obtained for the same 
system length L. In Fig.[n]we show nq{k, t) for a few fixed 
fc as a function of t. Almost independent of the system 
parameters, the temperature, and the considered mo- 
mentum the asymptotic behavior analytically described 
above (for L ^ oo) sets in very quickly and holds for 
times t < tr (oscillation frequencies and l/Vi decay). 
Just as for the energies the asymptotic values of nq{k,t) 
for large t and fixed k can be determined very accurately 
by averaging over an appropriate time interval. Similar 
to the behavior found for the energies the oscillations of 
nq{k,t) for fixed k become rather erratic if one exceeds 
ir (not shown here), but time averaging still gives very 
good agreement with the GGE result. Doubling the sys- 
tem size for the results shown in Figs. [S] and [5] does not 
lead to any changes on the scale of the plots and the 
curves can considered to be in the thermodynamic limit 
(as long as t < which is clearly the case in the figures; 
(«Fgc)ir ~ 6284). 

We can conclude that for the quench and the observ- 
ables of interest in the present work, which relax to a sta- 
tionary value, the latter is equal to the GGE prediction. 
The GGE expectation value differs from a thermal one 
(see the discussion in Sect. llVi in particular Fig. [2|). In 
the long-time limit the expectation values of the (subsys- 
tem) energies and the momentum distribution function of 
the fermions at every (fixed) k oscillate around the GGE 
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FIG. 6: (Color online) The same as in Fig.[S] but as a function 



of t for different k. For the momenta fc < 0, 1 



shown, 
arrows. 



The GGE expectation values are indicated by the 



result with a frequency AA((7min) and an amplitude de- 
caying as 1/Vt. In this respect the dynamics is rather 
simple — in particular compared to the one resulting from 
the fco-excitation discussed in the next section. 

We briefly comment on the relation of our calculations 
to those performed in the bosonization approach to the 
purely electronic Tomonaga-Luttinger model (in and out 
of equilibrium). For this model g-sums (-integrals) of the 
type appearing in Eq. (ISTt are usually performed analyt- 
ically after Taylor expanding the renormalized bosonic 
dispersion (here AA((7)) to linear order in q and assum- 
ing a particular g-dependence of the (= 1 — c^). It is 
generally believed that these steps do not alter the low- 
energy physics of the model. However, in Ref. l40l it was 
shown that this is not correct and the typical Tomonaga- 
Luttinger liquid exponents are affected by these approx- 
imations. In particular, this implies that the validity of 
the very interesting study of the relaxation dynamics of 
Ref. i is more restricted than it is realized by the au- 
thor(s). For the present model it is obvious that similar 
approximations would strongly alter the physics as the 
nonlinearity of AX{q) lies at the heart of our results. 



VI. THE DYNAMICS OF THE fco-EXCITATION 



Using the methods introduced above for the fco- 
excitation the expectation values of the observables stud- 
ied here have been computed in Ref. [1^ We emphasize 
that in this paper only the short-time behavior was inves- 
tigated while we are (mainly) interested in the long-time 
asymptotics. For completeness we here present the rele- 
vant expressions taken from Ref. ^251. 



given by 



cp 



p'^o (t) 



p'=o (t) 



p''0(t) 



= E«F9^ [4 + 4 + "^44 

q>0 ^ 

xcos{AA(g)i}] e(fco-9) , 

9>0 ^ 

X [1 - cos{AX{q)t}] e(fco - q) , (53) 
2n 



C/7 Sf) [ ^q} 



q>0 

X [1 - cos{AX{q)t}] e(fco - q) 



Lq ' ^ ' 



After the thermodynamic limit has been performed in the 
long-time limit the energy expectation values approach 
the GGE ones Eq. in an oscillatory fashion with the 
dominant frequency AA(qniin) and an amplitude decaying 
as 1/ \/t. This can be shown analytically applying the 
same methods as used for the quench. A typical example 
for the time evolution of the energies is given in Fig. 1 of 
Ref. m. 

For the phonon dynamics one obtains 



Nk, {q, i) = 2^ c2 s2 [1 - cos {AX{q)t}] e{ko - q) . (54) 



For the reason discussed in the last section the phonon 
momentum distribution function does not become sta- 
tionary. A plot of Nkg {q, t) similar to our Pig. |3] (obtained 
for the quench) is shown in Fig. 2 of Ref. [23. 

To compute {k, t) for the initial state containing 
a highly excited fermion in addition to the filled Fermi 
sea one has to bosonize four fermion fields instead of two. 
This leads to a rather involved expression for the fermion 
momentum distribution function as discussed in detail in 
the Appendix of Ref. [2^ In contrast to the quench this 
expression does not provide much analytical insight and 
in particular does not allow to analytically conclude, that 
in the long-time limit Uka {k, t) becomes equal to the GGE 
expectation value (at least after taking L — >■ 00 first). We 
thus have to rely on numerical comparisons. Therefore 
here we only give the equations needed for an iterative 
numerical calculation of {k, t) for finite 



no no 



-r min(r+s,r— ra) 



nko{kn,t) 



E E E 4:1,(0 

r=max(n,0) s=0 



1=0 



{ai-^Ht)yb["'\t)[bi-^\t)y , (55) 



for no > ri > and Uc = . The time dependent 
coefficients a\i^^ and are determined for m > 1, ^ G 



The time-dependence of the (subsystem) energies is Nq, and i = 0, m — 1 by 
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I 



[~-»m(t)/TO]^ ^(m-1) 



(56) 



J=0 



with the starting values (m = 1) 



k=Q 



k\ 



(1),,, _ [-Mt)]' [-u^m'-' , (57) 



(/-I)! 



-e{t) , Z e N 



6W(t) = 1 



The functions ui{t) and e{t) are defined as 



e{t) 



^—■iv-pt2TT/ L 



(58) 



In contrast to all our earlier expressions for the time de- 
pendence of expectation values both eigenenergies \±{q) 
enter explicitely and not only their difference AA(g). 
This leads to very rich dynamics. 

Figures [7] and [8] show two examples of nko{k,t) for 
a weak (F = 0.001) and an intermediate (F = 0.01) 
electron-phonon coupling at finite but large L. Only the 
part < k < ko is shown; for the time-evolution of ini- 
tially filled momenta see Figs. [10] and [11] At < = only 
the momenta fc < and feg are occupied (with weight 1) 
and no phonons are present. The "hot" electron decays 
into states with lower energy by producing phonons. The 
energy conservation in individual scattering processes is 
not sharp on short time scales which leads to broadened 
replicas of the initial /co-excitation at momenta fco — nq^, 
with n e N and — wo/wf- This is reminiscent to the 
derivation of Fermi's Golden Rule. In time-dependent 
perturbation theory the transition probability between 
two energy eigenstates with energy difference Auj is pro- 
portional to [sin(Awt/2)/(Aw/2)]^ which only at large 
t becomes an energy conserving (5-function. Thus the 
replicas sharpen for larger times. For our model they 
never reach the width of the original excitation before 
they get depleted again. Details of this dynamics are 
discussed in Ref. Two of such replica are visible in 
Fig. [7] for small g, while only a very broad feature ap- 
pears for intermediate couplings (see Fig. [5]). As soon 
as phonons are generated they couple to fermions in the 
filled Fermi sea (at A: < 0; not shown in the figures) and 
excite them to higher energies. This leads to a step like 
feature of width qb at small fc > and t > 0. For increas- 
ing time the sharp initial jump at A: = from 1 to is 
smoothened. In particular, the occupancies of the levels 
at A; ^ which are at t = will increase significantly. 
In Figs. [7] and [5] this is reflected in the continuing over- 
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FIG. 7: (Color online) The momentum distribution function 
of the fermions nkg{k,t). The parameters are F = 0.001, 



0.1, f = lO"'', and ko/qc 



0.5. The dark feature 



at fc = fco is the occupation of the intially filled level which 
strongly exceeds the scale of the z-axis. Note the scale of the 
«-axis. 
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FIG. 8: (Color online) The same as in Fig.[7]but for F = 0.01. 



all increase of nko{k,t) for small fixed k on time scales 
{{vFqc)t « 10^) at which saturation is already clearly 
established for the quench dynamics (see Figs. [S] and [5]). 
The initially filled level at fep is depleted but subsequently 
refilled on a time scale which obviously depends on the 
electron-phonon coupling (compare Figs. [7] and [S]). For 
large times (fco, t) approaches a small value in an oscil- 
latory fashion. For the quench the frequency with which 
nq(fc,i) oscillates (at fixed k) is independent of k and 
given by AA(gmin)- In contrast from Figs. [7] and [8] it is 
obvious that for the present initial nonequilibrium state 
(and on the same time scales as for the quench) different 
frequencies appear (compare the behavior at small fc and 
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FIG. 9: (Color online) The momentum distribution function 
of the fermions nkg{k,t) as a function of A; > for for fixed 
{vFqc)t = 10^ and different u = ix 10"^, 2 x 10"^, 10"^. The 
other parameters are F — 0.01, Q = 0.1, and ko/qc = 0.5. 
The inset shows nkQ{k,t) / v. 



k close to ko] see also Fig. [TT]). We conclude that the 
dynamics generated by the /co-excitation is significantly 
richer then the one found after a quench. 

Before comparing numerical data for nko{k,t) to the 
GGE result riGGE(fc) we have to gain a detailed under- 
standing of the subtleties of the thermodynamic limit 
for the present nonequilibrium initial state. We already 
noted that regardless of the electron-phonon coupling in 
the thermodynamic limit nGGE(^) is given by the ground 
state expectation value; the noninteracting step function. 
This result is consistent with our finding that the excess 
energies (both the GGE ones and the time-evolved ones) 
are not extensive. We only add the i-independent energy 
wpfco to the energy of the filled Fermi sea (the ground 
state). As the energy is conserved the three terms Eq. 
([5^ add up to vpko for all t > 0. Based on these consid- 
erations and our previous results we expect that 



lim lim nfc„ (fc, t) 

t—>oc L— )-oo 



e(fc) 



in accordance with the GGE prediction. This is con- 
firmed by Fig. [S] which shows nko{k,t) as a function of 
fc > for fixed (uf9c)^ — 10^ and different = 4 x 10~^, 
2 X 10~^, 10~^. In the inset we plot nk„{k,t) /v and the 
collapse of the three curves indicates that {k, t) for 
fc > vanishes as 1/L. The same holds for 1 — nko{k,t) 
at A: < 0. 

After clarifying the behavior in the thermodynamic 
limit we now compare the time evolved and GGE dis- 
tribution functions for fixed L < oo. In Fig. [TU]we show 
the momentum distribution function as a function of k 
for a few fixed t, the same model parameters as in Fig. 
[5] (quench), and ko/qc — 0.5. Up to very small momenta 
|fc| ^ Qc (see the inset) we find convergence towards the 
GGE prediction for times t < t^. As argued above the 
changes during the time evolution for the small momenta 
are large (from a step function at t = to a continuous 
function at large t) and convergence cannot be found for 
times t <t^. Note that because of the slow convergence 
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FIG. 10: (Color online) Momentum distribution function of 
the fermions nkg(k,t) as a function of k for different t. The 
GGE distribution obtained for the same system size is shown 
for comparison. The parameters are F — 0.03, f2 = 0.1, v = 
10"^ and ko/qc = 0.5. 
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FIG. 11: (Color online) The same as in Fig. 1101 but as a 
function of t for different k. For the momenta fe < 0, 1 — 
rifcQ {k, t) is shown. The GGE expectation values are indicated 
by the arrows. 



the times shown in Fig. [10] are much larger than those of 
Fig. [5] even though the model parameters are the same 
in both figures. Finally, in Fig. [Tl]we show nko{k, t) for 
a few fixed A: as a function of t. The GGE expectation 
values are indicated by the arrows. For k sufficiently far 
away from the initial step at fc = it is again evident that 
nko{k, t) approaches ?t.gge(A). In that sense a "quasi sta- 
tionary state" (finite L) is reached with a time averaged 
expectation value (over "large" times smaller than t,.) 
which can accurately be described by the GGE predic- 
tion and which is different from the step function reached 
in the thermodynamic limit. In complete analogy to our 
earlier findings for the quench dynamics extending the 
time beyond tr leads to enhanced oscillations with a time 
average which is still close to nGGE(A;). 

We note in passing that for the strongest possible 
electron-phonon coupling allowed by stability F = f2, 
A_(q) = for all q < gc- A similar situation is discussed 
in Ref. H. In this case the amplitude of the oscillation in 
nko{k, t) for fixed k does not decay. The energies still ap- 
proach a constant large time limit but with a 1/t instead 
of a l/\/i decay 4i 
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We did not succeed in extracting analytical results for 
the long-time asymptotics from the rather involved finite 
L expression for nko{k,t) presented in the Appendix of 
Ref. I25I . A numerical analysis for curves which at suffi- 
ciently large times (but still t < t^) oscillate around the 
finite L GGE prediction (e.g. the k/qc~ —0.4, —0.2, and 
0.3 curves in Fig. [TT|) shows that the amplitude decays 
faster than in the case of a quench, that is faster than 
l/\/t. This has to be contrasted to the observation that 
it takes much longer times in case of the fcg-excitation 
before any asymptotic behavior sets in (compare Figs. [S] 
and llip . From the numerical data it turned out to be im- 
possible to extract an analytical form of the amplitude 
decay (e.g. power law with an exponent x larger 
than 1/2). As already noted when discussing Figs. [7] and 
|5] the oscillation frequency of riko (^1 '"^^ fixed k seems 
to dependent on k. This becomes even more evident in 
Fig. [TTJ Besides this the oscillation frequency depends 
on the model parameters F and fi. We did not succeed 
in extracting a clear picture for the dependence of the 
frequencies on these three parameters. In addition, it 
is not obvious how to relate the numerically determined 
frequencies appearing in nko{k^t) at fixed k to the mode 
energies X+(q) and A_(g). This again exemplifies that 
the time evolution of the fermionic momentum distribu- 
tion function for the fco-excitation is rather involved and 
much richer than the dynamics encountered above. 



state dependent generalized Gibbs ensemble. We found 
that for all observables of interest to us the asymptotic 
long-time limit of expectation values which become sta- 
tionary are equal to those obtained from the appropriate 
natural GGE. We discussed that the momentum distri- 
bution functions of the GGE differ from those obtained 
within the canonical ensemble. While the (excess) en- 
ergies and the electron momentum distribution function 
become stationary, the phonon momentum distribution 
at fixed momentum oscillates with a constant amplitude 
even at large times. This can be understood from the lin- 
ear relation between the eigenmode ladder operators and 
the phononic ones. Long-time convergence of the ener- 
gies and the fermionic momentum distribution function 
in the strict sense is only achieved if the thermodynamic 
limit is taken first. Going beyond this we showed that 
the GGE predictions for these expectation values agree 
to time averaged values taken at large times (even larger 
than the characteristic scale tr — L/v-p). For most of 
the studied situations we were able to analytically de- 
scribe how the GGE prediction is reached. We found 
that the asymptotic limit is only reached following a 
power law (instead of exponentially). The dynamics of 
the fermionic momentum distribution function resulting 
from the /cQ-excitation turned out to be much richer than 
the one found for other observables and for the interac- 
tion quench. 



VII. SUMMARY 

We studied the relaxation dynamics of an exactly solv- 
able electron-phonon model out of two initial nonequi- 
librium states. This way we added a continuum model 
containing both fermions and bosons to the list of quan- 
tum many-body models recently investigated and con- 
sidered other nonequilibrium situations than the param- 
eter quench. Our model features a natural set of con- 
stants of motion, with as many elements as degrees of 
freedom. This set can be used to construct the initial 
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